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Abstract 
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using chiral perturbation theory to next-to-leading order {0{e'^Q)) at energies, w < m^^, where m^r 
is the pion mass. Our results show that the unpolarized differential cross section can be used to 
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are sensitive to different linear combinations of the four neutron spin polarizabilities. 
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I. INTRODUCTION 



The response of an object - that has sub-structure - to a quasi-static electromagnetic 
(EM) field is characterized in terms of quantities called polarizabilities. For example, when 
an object, having charged constituents, is placed in an electrostatic field the centers of 
positive and negative charge separate resulting in an induced electric dipole moment, p = 
aE, where E is the external electrostatic field. The strength of this response, a, is defined as 
the "electric polarizability" . Similarly, we can imagine a magnetic dipole moment /I = f3H 
being induced in the presence of an external magnetic field, H, and this magnetic response is 
quantified in terms of the magnetic polarizability, (3. If an object has intrinsic spin, then in 
the presence of an EM field the orientation of the spin may be altered. This is the spin-dipole 
response; in fact, spin responses are quantified in terms of four spin polarizabilities. Here, in 
this work we shall focus on the electric, magnetic and the spin polarizabilities of the neutron. 
Investigating the neutron polarizabilities would contribute immensely to our understanding 
of the neutron structure and also proton- neutron as isospin-doublet. We investigate neutron 
polarizabilities by calculating elastic Compton scattering off '^He because there are no free 
neutron targets and ^He behaves as an 'effective' neutron (as we shall see later). 

For the following discussion we consider both the proton and the neutron, which we 
commonly call nucleon. Using the concept of an effective theory, the most important contri- 
butions for the interaction of a nucleon with an EM field can be identified. The interaction 
Hamiltonian can then be formulated in terms of the EM fields and nucleon EM moments. 
The symmetries of the EM interaction ensure that only those terms need to be considered 
that obey the following rules [l|: 

1. Only terms quadratic in A (the EM vector potential) are allowed. 

2. The Hamiltonian is gauge invariant, 

3. a rotational scalar, and, 

4. P (parity), T (time-reversal) even. 

Defining uj as the photon energy, the leading-order response in uj for a nucleon in an EM 
field is given by- 

Hi^ = + e<^, (1) 

eff 2M ^ ^ 

where A and are the EM vector and scalar potentials. At the next order the following 
terms appear in the Hamiltonian (see e.g. [H, Uj): 



eff 2 



aEiE'^ + PMiB^ ■ (2) 



Here, oei and (3mi are the electric and magnetic dipole polarizabilities and are the most 
prominent of the nucleon polarizabilities. This Hamiltonian (Eq. ([2])), however, includes only 
terms of second order in u and can be extended, as long as the aforementioned conditions 
are respected. Apart from aEi and Pmi there exist other polarizabilities, such as the spin- 
dependent polarizabilities. Therefore, we can modify and rewrite Eq. ([2]), extended to the 
spin dipole polarizabilities [H, [2|: 



aEi E^ + Pmi 5 + lEiEi a- E X E + 'Jmimi a- B x B 
2 7mi£;2 CTi Eij Bj + 2 ^eim2 cTi Bij Ej] , (3) 



2 



with a denoting the intrinsic nucleon spin. The various subscripts in a, /3 and 7 denote 
different multipoles of the nucleon's response to the external EM field and 



(4) 
(5) 



In this work, we use 71, 72, 73 and 74 [3| to represent the four spin polarizabilities and these 
are related to the 7's in Eq. ([3]) as - 



lElEl 
IMIMI 
7M1E2 
1E1M2 



-71 - 73, 

74, 

72 + 74, 

73- 



(6) 



An electromagnetic probe by nature, Compton scattering captures information about the 
response of the charge and current distributions inside a nucleon, and hence the polarizabil- 
ities, to a quasi-static electromagnetic field. To lowest order in photon energy (C(cc;°)), the 
spin-averaged amplitude for Compton scattering on the nucleon is given by the Thomson 
term 

Amp = —e ■ e', (7) 



M 



where Ze^ M represent the nucleon charge and mass respectively and e, e' specify the po- 
larization vectors of the initial and final photons respectively. At the next order in photon 
energy, contributions arise from electric and magnetic polarizabilities — ue and (3m — which 
measure the response of the nucleon to the application of quasi-static electric and magnetic 
fields. The spin-averaged amplitude is expressed in the laboratory frame as: 



Amp 



e ■ e 



2„2 



+ ujuj' AnaE ] + e x k ■ e' x k' (jJuAttPm + 0{uj^ 



(8) 



Here, k^ = {uj,k), k'^ = {uj',k') specify the four-momenta of the initial and final photons 
respectively and k and k ' are unit vectors associated with the photon momenta. When we 
square the amplitude ([8]), the associated differential scattering cross section on the proton 
is given by 



da 

In 



( 



, -V 

V47rM J 
AttMulu' 



{oe + /5a/) (1 + cos( 



f + ^{aE-PM){l-cose)^ 



(9) 



From Eq. ([9]) it is evident that the differential cross-section is sensitive to {ue + Pm) at 
forward angles and to (a^ — Pm) at backward angles. 

It is essential to mention here that for the purpose of this work, oe and Pm are the so- 
called Compton polarizabilities and measure the "true deformation" effect on the nucleon. 
In other words, Compton polarizabilities measure the distortion effects of charge and current 
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distributions inside a nucleon subjected to an EM field. In the non-relativistic limit, these 
reduce to the static electromagnetic polarizabilities described in Eq. ([2]). For example, the 
static electric polarizability, a£;i = Q;£; + f^|3], where k is the anomalous magnetic moment 
of the nucleon and M is its mass. From here on we shall drop the subscripts E and M as well 
as the bars over a and (3 from the Compton polarizabilities. Therefore all polarizabilities 
discussed from now on are the Compton polarizabilities. At the next order in uo beyond that 
considered in Eq. ([H]), one can access the spin polarizabilities (7's) of the nucleon. 

The last fifteen years or so has been witness to a concerted effort on the part of both 
theorists and experimentalists to comprehend the structure of nucleons as manifested in 
Compton scattering. Through various experiments on the proton, its electric and magnetic 
polarizabilities have been effectively pinned down. There has been an avalanche of experi- 
mental data on unpolarized Compton scattering on the proton at photon energies below 200 
MeV js], IgI, 0, IsI, @]. The current Particle Data Group (PDG) values for the proton are: 

Up = (12.0 ±0.6) X 10-^fm^ 

Pp = (1.9 ±0.5) X lO^^fml (10) 

Most of the extractions of the proton electromagnetic polarizabilities were performed via 
a dispersion relation approach by a multipole analysis of the photoab sorp tion ampli- 
tudes [lo|, 11]. Calculations have also been done up to 0{e^Q^) in xPT (ij [H, Q 

and 



these calculations give a very impressive description of the data for u, ^/\t\ <200 MeV. 
Note that it is necessary to measure only one of a and (3. In practice, one can then use 



the dispersion sum rule (Baldin sum rule) derived from the optical theorem [15 



<*(cu) 



01 + 13 = —: ^^^dcu, (11) 



that relates the sum of the polarizabilities to the total photoabsorption cross-section in order 
to extract the other. In Eq. (ITT!) . cr*°*(ci;) is the total photoabsorption cross-section and Uth 
is the pion-production threshold. For the proton 



Op 4- /3p = (13.8 ± 0.4) X 10"^ fm^ (12) 

There are several evaluations of the dispersion relation in Eq. (ITT]) for the neutron [2I, 
13]. However, there is discrepancy between the numbers extracted. A recent extraction of 
the sum rule by Levchuk and L'vov [13] using relevant neutron photo-production multipoles 
gives 

a„ + /5„ = (15.2±0.5) X lO^^fml (13) 

Because of the discrepancy in the extraction of the neutron sum rule, as far as the neutron 
polarizabilities go, we still strive to extract precise numbers for the polarizabilities. Since 
neutrons are very short-lived, we lack free neutron targets and this poses a serious handicap 
in the study of neutron structure. The drawback caused by not having free neutron tar- 
gets encouraged the community to look at other avenues to extract information about the 
neutron. 

There were attempts to measure neutron polarizabilities via scattering neutrons on lead 
to access the Coulomb field of the target and examining the cross-section as a function of 
energy. Currently there is much controversy over what this technique gives for a„. Two 
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experiments using this same technique obtained different results for a„ (Refs. |18l. Il9|). 



ar 



(12.6 ± 1.5 ±2.0) X 10"^ fm^; 
(0.6 ±5.0) X 10"^fml 



(14) 
(15) 



Enik et al 20j revisited Ref. [l8| (Eq. (IT^ ) and recommended a value of a„ between 7 x 10 ^ 
fm^ and 19 x 10~^ fm^. 



In addition, quasi-free Compton scattering from the deuteron was measured at SAL [21 
for incident photon energies of = (236 — 260) MeV and one-sigma constraints on the 
polarizabilities were reported to be: 



(7.6-14.0) X 10"^ fm^ 
(1.2 - 7.6) X 10"^fml 



(16) 



Another quasi-free experiment was performed at Mainz 2^ and the difference in the polar- 
izabilities was obtained: 



a, 



-Pn = (9.8 ± 3.6(stat) ± 2.2(model)+^:}(sys)) x 10"^ fm 



(17) 



Until now, however, there have been only limited efforts to measure nucleon spin po- 
larizabilities. The only ones that have been extracted from the experimental data are the 
forward and backward spin polarizabilities. The backward spin polarizability is defined as: 



Ttt = 7i ± 72 ± 274. 
The neutron backward spin polarizability was determined to be 

j^n = (58.6 ±4.0) X 10-^fm^ 



(18) 



(19) 



from quasi- free Compton scattering on the deuteron [22|. This experiment used the Mainz 48 
cm X 64 cm NaT detector and the Gottingen recoil detector SENECA in coincidence. For 
comparison, an experiment 2^ using the same detector set-up extracted a proton backward 
polarizability value ranging from (-36.5 to -39.1)xl0~^ fm*^. These values are consistent 
with earlier Mainz measurements p, 24 , 25l . 



The theoretical prediction for 7,^^ from xPT is —36.7 x 10~^ fm^ 26|, |27| and these 
calculations were done up to the next-to-leading order with the A-isobar as an explicit 
degree of freedom. The prediction xPT for 7^^ is 57.4 x 10~^ fm^ |27i]. Both predictions are 
in agreement with the experimental numbers to within their quoted uncertainties. 

The forward spin polarizability, 70, is related to energy- weighted integrals of the differ- 
ence in the helicity-dependent photoreaction cross-sections {cri/2 — 173/2). Using the optical 
theorem one can derive the following sum rule for the forward spin polarizability [isl. [28|: 



7o = 7i - (72 ± 274) 



47r2 



0"l/2 — 0"3/2 



da;, 



(20) 



where Uth is the pion-production threshold. The following results on 70 were estimated using 
the VPI-FA93 multipole analysis [2^ to calculate the integral on the RHS of Eq. fl20l) : 



7op 

70n 



-1.34 X lO^Mm 



-0.38 X 10"^ fm 



(21) 
(22) 
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The xPT prediction for 7op is —2.1 x 10~^ fm^ 2^ and for 7on is consistent with [27| . 

Let us now summarize the "state of the nucleon polarizabihties" . The quantities ap and 
f3p are well known, an and f3n are also known. However, looking at the numbers in Eqs. (fT3|l - 
( TT7|) . it is evident that the extractions of a„ and /3„ have large error bars. Strikingly, there 
is no consensus on the numbers for 7ip . . . 74^ or on 7i„ . . . 74^. Of these, jip . . . 'j^p are not 
known due to a lack of experimental data. But it is clear that the extraction of a„, and 
7in • • • 74n require efforts from both the theoretical and the experimental community. 

As there are no free neutron targets, the theoretical community has directed focus on 
elastic Compton scattering from light nuclei to access the neutron polarizabihties. The 
main advantage of such a reaction process compared to the more complex n-Pb scattering 
is that the theoretical analysis of the few-body reactions is much better controlled. The 
lightest nucleus is the deuteron and elastic •yd scattering has been studied with the intent of 
extracting information about the neutron polarizabihties from unpolarized and polarization 
observables. It should be noted, however, that processes with A > 1 pose a different kind of 
challenge because they require an understanding of the inter-nucleon interaction. Deuteron 
structure is governed by the NN interaction and hence, when analyzing 7d data, one has to 
include the effects of photons coupling to mesons being exchanged between nucleons ('two- 
body currents') over and above the single-nucleon 7N amplitude. There exist calculations 



for 7d scattering in the framework of conventional potential models [17|, l30|, l31|, l32, l33|, [34 
All of these calculations use a realistic NN interaction and usually give a good description of 
the data. However, the predictions of these calculations depend somewhat on the NN model 
chosen. xPT on the other hand allows a model independent framework to systematically 
build the theory for 7d scattering. The amplitude for coherent Compton scattering on the 
deuteron has been calculated to C(e^Q^) by Beane et. al. [H, At this order there 
are four new parameters in the theory (two of which can be fit to the jp data) and since 
the deuteron is an isoscalar, the nucleon isoscalar polarizabihties were extracted from the 
above experimental data by fitting the other two parameters. The isoscalar polarizabihties 
are linear combinations of the proton and neutron polarizabihties, = '■"p^""-' and /3n = 
i&^i^. Beane etal. obtained- 

ttTV = (8.9 ± l.5)tol X 10"^fm^ 

(3n = (2.2 ± 1.5)+^:^ X 10"^ fm^ (23) 

Subsequently, next-to-leading order calculations for deuteron Compton scattering have been 
performed that include the A as an explicit degree of freedom and also incorporate the 
low-energy NN rescattering contributions in the intermediate state [H, [s^, 33|. Thus, so- 



phisticated calculations on 7d scattering already exist and with more experimental data 
soon to come from MAXLab in Sweden there is hope that and /?„ can be extracted with 
good precision from that experiment. 

Recent experimental advances (for example technology to polarize targets) have made it 



possible to venture into the area of polarization observables. Hildebrandt etal. |37|, |38| have 
also studied polarization observables for Compton scattering on the nucleon (both proton 
and neutron) with a focus on the spin polarizabihties. Their calculation involves free nucleons 
and may be beneficial for understanding processes that involve quasi-free kinematics like 
'jd — >• ■ynp. Choudhury and Phillips i40] focused on ^d jd and jd 'yd to analyze the 
effects of the electric, magnetic and the spin polarizabihties of the neutron. They reported 
that one of the double-polarization observables was sensitive to a linear combination of the 
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neutron spin polarizabilities and the deuteron Compton scattering observables alone were 
not enough to extract the neutron polarizabihties. Calculations for polarization observables 
for deuteron Compton scattering with explicit A-isobar and low-energy resummation have 



also been completed and the results are forthcomi ng 14 1 



These studies [11, Q [13, H Hll, S H S, SrS, H S, M, 0, El, 113 make it 
evident that not all of the neutron polarizabilities can be extracted from deuteron Compton 
scattering. This means that several different observables or combinations of observables 
from different reactions are necessary to effectively pin down the neutron polarizabilities. 
To this end, physicists should focus on other alternatives, for instance, quasi-free kinematics 
in deuteron Compton scattering, or other nuclear targets to extract information about the 
neutron. Such parallel calculations to study the effect of the same neutron polarizabilities 
combined with experiments can build confidence in the extraction of the polarizabilities. In 
principle, one should extract the same numbers, no matter what the process or the target is. 
Thus, these parallel studies would reassure us about our understanding of "nuclear" effects. 

One alternative target is the polarized ^He nucleus. This nucleus has the nice property 
that the two proton spins are anti-aligned for the most dominant part of the ^He wavefunc- 
tion, the principal 's-state'. Since approximately 90% of the wavefunction is given by this 
state, the spin of ^He nucleus is mostly carried by the unpaired neutron alone. Fig. [1] shows 
the configurations for the principal 's-state' and for one other possible contributions to the 



He wavefunction [43|, . This paper reports our calculations for Compton scattering on 



a '^He nucleus. Note that, for 7^He scattering we now have to understand the interplay of 



P P n 



P P n 




Principal s-state: 90% 



d-state: 8% 



FIG. 1: The polarized ^He nucleus. 90% of the time the nucleus is in the principal s-state 
configuration. 



the three nucleons over and above the two-nucleon effects. 

In this work, we calculate elastic Compton scattering on ^He to 0{e^Q) and focus on 
specific observables so as to construct road-maps to extract the neutron polarizabilities, 
a„, Pn and 7i„ . . .74„. It should be mentioned here that these are the first ^He Compton 
scattering calculations in any framework. These calculations include the processes 7^He^ 
7^He and 'j^He — > 7^He and are necessarily exploratory. They are expected to provide 
benchmark results for elastic Compton scattering on ^He. As such, our results will serve 
to attract further explorations of these observables at Oje^Q"^) and beyond. The main 



conclusions of this work have already been reported in Ref. 




This paper is intended 
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to be a comprehensive discussion on the calculations that lead to these conclusions. 

This paper is organized in the following manner. Sec. [TTl lays out the structure of our 
calculation. In Sec. [TTllwe describe the various observable that we focus on and then Sec. |V] 
we report our results. In order to make this work comprehensive, a discussion on inherent 
sources of uncertainties that may affect the final result is provided in Sec. IVIII Finally, 
in Sec. IVIIIt the calculation is summarized and also future directions for calculations of 
7'^He scattering are pointed out. 



II. ANATOMY OF THE CALCULATION 

The irreducible amplitudes for the elastic scattering of real photons from the NNN system 
are first ordered and calculated in HB^PT . Then these amplitudes are sandwiched between 
the nuclear wavefunctions to finally obtain the scattering amplitudes (see Fig. [2]). This 
amplitude can be written as- 

M = (24) 

Here, |\E'j) or |\I'/) are the ^He wavefunctions that are anti-symmetrized in order to take 
into account that the nucleons are identical fermions. We will employ ^He wavefunctions 
obtained from Faddeev calculations in momentum space. We first calculate a specific Fad- 
deev component for ^He which is related to the fully anti-symmetrized wavefunction 
I*) H by 

|vI/) = (l + P)|^)3 (25) 

where P = -P31-P12 + -P32-P12 is the sum of the cyclic and anti-cyclic permutation operators and 
the numbered subscripts denote the nucleon number. In our case, since we are calculating 




FIG. 2: The anatomy of the calculation. The irreducible amplitude is computed in HBxPT and 
sandwiched between external ^He wavefunctions to give the matrix element for Compton scattering 
on ^He . 



only up to 0{e^Q) in HBxPT , the operator O consists of a one-body part 



O 



IB 



0^^(l)+0^^(2)+0^^(3) 



(26) 



and a two-body part 



O 



2B 



02^(1, 2) +0^^(2,3) + 0^^(3,1). 



\2B 



2B/ 



(27) 



The mechanisms that contribute to the two different pieces are shown in Figs. [3l H] and O 
The structure of the 7N and the 7NN amplitudes follows Refs. [l^, 42] and they are shown 
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in Apps. 1X1 and [Bl In the notation used in Eq. 0261] . O^^(a), a = 1 ... 3 represents the one- 
body current where the external photon interacts with nucleon 'a' and in Eq. (^7^, 0'^^{a, b) 
represents a two-body current where the external photon interacts with the nucleon pair 
'(a, 6)'. In short, the different terms in Eqs. fl2Bl) and fl7r|) represent different permutations 
of the three nucleons inside ^He . Note that, we do not have any three-body currents because 
they appear only at 0{e^Q^) in the xPT power-counting. 




FIG. 3: Tree level diagrams for the one-body amplitude. 



If we now take into account the identity of the three nucleons, we can rewrite Eq. 

as- 

>1 = 3(%|6i^(l) + 62^(1,2)1^,) (28) 



For our calculations we draw upon the work of Kotlyar et al j47|. They developed an ap 



proach for calculating matrix elements of meson-exchange current operators between three- 
nucleon basis states. They applied this approach in studying ^He photo- and electro- 
disintegration. The three-nucleon basis states were expressed in a jj-coupling scheme (final 
configuration) and a three-nucleon bound state (initial configuration). These matrix ele- 
ments were then expressed in terms of multiple integrals in momentum space. Since they 
studied photo- and electro-disintegration, they had only one external momentum transfer 
vector q that they could choose to align in a preferred direction in order to simplify the 
calculations. Here in contrast, since we have an incoming and outgoing photon with differ- 
ent three-momenta, we do not have the same liberty and hence we needed to extend their 
prescription. Also, we choose to calculate the one-body currents too in the two-nucleon 
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✓ V 

/ \ 
/ \ 

< ^=^ 

FIG. 4: Contribution to the one-body amplitude at 0{e'^Q). These diagrams contain one pion 
loop. 

spin-isospin basis and hence it is more convenient to express Eq. ( |28|1 as- 

A< = 3(v[/^|^(di^(l) + 0^^(2)) +02^(1,2)1x1/,) =3(vl/^|6(l,2)|xl/,). (29) 

The advantage of this formulation is that the structure of the calculation is similar for both 
the one-body and the two-body parts. The next step is then to actually calculate the matrix 
elements. To do this, we project the state-vectors and the operators on to the basis- 

bl2P3tt) = bl2P3)|(^12Sl2)jl2(/3^)j3(j3jl2)^Mj)|(ti2^)TMT)3 = 1^12^3) l^j) |aT)3- (30) 

Here, pi2 and denote the magnitude of the Jacobi momenta of the pair "(1,2)" and the 
spectator nucleon "3" and this choice ensures that the set of basis states are complete and 
orthonormalized. The subscript '3' represents the choice of the Jacobi momenta (Fig. [6]) 
defined as- 

Pi2 = h) 

P3 = li2h -h-h). (31) 
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FIG. 5: Contributing two-body diagrams at 0{e^Q). Permutations are not shown. 



i.e. nucleon '3' serves as the spectator. The total angular momentum of the nucleus J is a 




FIG. 6: Jacobi momenta for the three-nucleon system. 



result of coupling between ji2 and j'a, the total angular momenta of the '(1,2)' subsystem and 
the spectator nucleon '3' respectively. The orbital angular momentum I12 and the spin S12 
of the two-body subsystem '(1,2)' combine to give ji2 and similarly, ^3 and S3 = | combine 
to give js. ti2 is the isospin of the two-nucleon subsystem and it combines with the isospin 
of the spectator nucleon to give the total isospin T and Mt is simply the projection of T. 
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Since we are concerned with the ^He nucleus- 

\TMt) = I--) 
I 1/ I22/ 

\JMj) = \^Mj). (32) 

The basis states as defined above are not completely anti-symmetrized. The antisymmetry 
is carried by the complete wavefunction. However, for the two-body subsystem quantum 
numbers the antisymmetry within the subsystem leads to the additional constraint that 
h2 + S12 + ti2 = 2n + 1, (n = 0, 1, 2, . . .)• 

Before delving into the matrix element calculations, let us first consider the kinematics 
of Compton scattering from a ^He nucleus with momenta assigned in this way. 



Kinematics 



Since we are interested in two-body currents at the most, the kinematics of our problem 
becomes a little simplified. In our calculations we shall be working in the 7 ^He cm. frame. 
As we can see in Fig. [TJ nucleon "3" serves as a mere spectator to the 7NN —>■ 7NN scattering 
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FIG. 7: The kinematics of the calculation in 7^He cm. frame. The vectors k {—k) and k' {—k') 
are the three-momenta of the incoming and outgoing photon (^He) in the 7^IIe cm. frame. 

process and this means that in our calculations we can treat the "(1,2)" and "3" momentum 
spaces separately. Moreover, since the third nucleon does not partake in the interaction, we 
have a three-momentum conserving delta-function from- 



P3 - 3 = P3-3 



P3 



P3 + 



k 

3 

(k' - k) 



P3 + 



1 

3' 



(33) 
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FIG. 8: The co-ordinate system. 




FIG. 9: Compton scattering - before and after picture. 



where, q = {k' — k). Hence, in principle, any integral in p^, ' space can be eliminated by using 
this simplification. However, as is evident from Fig. [71 there are still three unknown three- 
momenta and hence we have to perform a nine- dimensional integral to calculate the matrix 
elements. We shall see later that the matrix-elements pertaining to the one-body currents 
are further reduced to six- dimensional integrals because of an additional three-momentum 
conserving delta-function. 

Before we go further into the discussion of the calculation, let us first mention that in 
the subsequent descriptions the following convention for the co-ordinate system is used (see 
Fig. [8]). The beam direction is defined as the z axis. The x — z plane is the scattering 
plane with the y axis being normal to it. Fig. [9] shows the before and after picture for 
Compton scattering from a target. The wiggly line depicts a photon, the red circle depicts 
the target. As mentioned earlier, the scattering takes place in the x — z plane with 6 defining 
the angle of scattering in the center of mass (cm.) frame. For linearly polarized photons, 
the polarization vectors in the initial state can be along the x or the y axis, i.e., e = x or y. 
For circularly polarized photons in the initial state, the beam helicity. A, can be ±1. The 
photons are right circularly polarized (RCP) if the beam helicity, A = +1 (e_|_ = — or 

left circularly polarized (LCP) if the beam helicity, A = — 1 (e_ = ^7^). 
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B. The Matrix Element 



Focusing back on our matrix element calculation, we now introduce complete sets of 
momentum states to obtain- 

MiM'j.Mj) = 3{^filM'j)\dil,2)\^,i^;Mj)) 

= 3 j d^p,2' j dW j d^p,2 j d'p, Yl 

o,a',mi2,m'j^2 

M'j) |pi2 % V) {pu % V|0(1, 2) \pi2P3a) {puPsal^ Mj)) 
d^Pu ' j d^P^ ' j d^pu j d^p3 

'P{Pi2\P3',a')^{Pi2,P3,a) 

a,a',mi2,»Tij2 

{Pl2'p3']aT\0{l,2)\pi2P3]aT)yh2,sr2jl2,mAPl^^^ (34) 

where |\E'j(|; Mj)) and (^^/(l; Mj)] denote that the '^He nucleus is in a specific isospin and 
total spin projected state. In Eq. flMl) . we define 

yi,s,j,m{p) = ^(/,m-m5,s,ms|/,s,j,m)Y;,„_„^(p)|sms) (35) 

ms 

and 

V5(Pi2,P3,a) = {ji2,mi2j3,Mj -mi2\ji2,j3,J,Mj) 

(ti2, mti2, - , mtgltia, - , -, -)(pi2P3tt|^)- (36) 

where the ^He wavefunction (pi2P3a|^) is independent of Mj. 

For a consistent calculation of the matrix element, the nuclear interaction used to generate 
the ^He wavefunction should be based on the same chiral effective field theory approach 
that we used for the derivation of the operators. Strictly speaking, the available nuclear 
interactions have been derived in slightly different frameworks. Therefore, in order to get an 
idea of the dependence of our results on the choice of the wavefunction, we have employed 
several nuclear interactions to obtain the ^He wavef unctions. These are the chiral NLO 



interactions of Re fs. 1481 . |49|, the chiral Idaho N^LO interaction |50j and the high precision 



model interaction [51[ . In order to correct for the underbinding of the He system when only 



NN interactions are used, we augment these interactions also by NNN interactions [52, l53 
as specified later. Based on this range of interactions, the binding energy of ^He is between 
6.89-7.83 MeV (the experimental one is 7.72 MeV). It turns out that a rather small number 
of partial waves is sufficient to sufficient to achieve convergence for the Compton scattering 
matrix elements. We found that the one-body (two-body) matrix elements are converged to 
within 0.1% (0.2%) using partial waves with ji2 <2 (ji2 <1)- 

In the sum in Eq. flM|) . we can use the fact that the isospin projection of the third particle 
remains unchanged, as does the total isospin projection, to introduce the kronecker-delta 
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Smt[2,mti2- Eq. is a 12-dimensional integral and it can be reduced to a 10- dimensional 
one by simply implementing the delta-functions over the angular parts of p^'. Since we are 
separating the angular and magnitude part of the three-momenta, from here on we shall 
drop the 'vector' sign in the momenta; for example {ps'l = Ps'- This leads to- 

MiM'j,Mj) = 3 Yl 

j Pi2'^dpi2' j Pl2dpi2 j dp'^ j pldp3 
(p{pi2',P3',a)^iPi2,P3,a) 

j dpu' j dpi2yp^^^s[^j:^^,m[,iPl2'){Pl2'; t[^mti2\d{l, 2)\pi2] ^12"^^12) 
3^/l2,.l2,ii2,mi2(Pl2) j "^P^yilf^M'j-m',, + |^ 

S(P3' - 1^3 + 3^,3,1 j3,M,-mi2 (^3) 

il2 >il2 1™! 2 1™12 iS'l2 >*12 i''l2 i'l2 J3 lis ,^3 ,^3 ,mtl2 
j Pl2''^dpi2' j Pi2rfpi2 j pldp3Lp{pi2,P3,a) 
^2(^12, P12 '-Jl2, l'l2^ S12, S'i2, jl2, jl2' "^12, ^2' "^^12) 



dpz V(Pi2 P3 a')l3{p3, P3; /a, /g, jg, jg, mi2, m'la, Mj, Mj). (37) 



The integral 



X2(pi2, Pl2 lu, l'l2, -312, S'i2, jl2, 7^2) "^12, "^12) "^^12) = 
y #12'! dpi2yp^^,s[,,r,,,rn'JPl2') 

{Pl2 t'i2"^^12|C'(l, 2) \pi2; tl2"^il2)3^/i2,si2 Ji2,mi2 (P12) (38) 

is a four- dimensional integral and is computed numerically. By using the delta-function we 
can reduce the dimensionality of 

2^3(^3, P3, I's, is' ^12, Mj, M'j) = 

/ "^P^^ihi^'M^j-m',, + 1)^{P^'- '^"S + ||) 3^,3, l,,3,M,-„^l2 (P3) (39) 

which is two-dimensional, thereby making the final integral a nine-dimensional one. We 
shall first talk about how I3{p3, p'3] I3, 13, jsyjs, 17112, '"^'12? Mj) is reduced and then focus 
on X2(pi2,Pi2 h2, ^12, ■S12, s'i2, ji2, ji2, f^u, 'm'i2i f^tu)- Howcvcr, even at that point our job is 
far from being done because the operators also have to be evaluated in this same basis. We 
shall discuss this procedure right after we finish discussing how the integrals were calculated. 
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1. Reduction of Integral I-i 



As noted in Eq. fl37j) . X3 is a two-dimensional integral involving the angular integration of 
the momentum P3. However, we can translate the delta-function in |p3'| to one that involves 
the azimuthal angle of using the following steps. 



2ps'S (p3'^-|P3 + |r 



12 



2 1 



2p3q 



(pPq) ■ 



(40) 



Since, we work in the co-ordinate system defined in Sec. Ill Al k = (0,0,1) and k' = 
(sin ^, 0, cos ^) and this means that q is in the x — z plane or (pq = 0. Then defining 
Ps = (sin 6*3 cos 03, sin 6'3 sin 03, cos 6^3) we can rewrite Eq. fHU]) as- 



3P3' 



p^q sin ^3 sin 9q 
Sps' 1 



5 I cos 03 

1 



P3q sin 03 sin 9q I sin , 



Z3 — cos 6*3 COS 9q 

sin 6*3 sin 9q 



5(03 - 03) + 5(03 + 03) 



(41) 



where, 



Using Eq. 
5 ( P3 ' - \P3 + 



Z3 



cos 03 



P3 



12 



pI 



Q_ 

9 



2p3g 
3 



Z3 — COS 6*3 COS 9q 

sin 9s sin 9q 



(42) 



we can now eliminate the 03 integral in X3. The delta-function 
further sets the bounds of integration for both the P3' and the ^3 in- 
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tegrals. Implementing the above steps, we can now express the last line of Eq. (1371) as- 
dps V(Pi2 P3 a')2^3(P3, Ps, h, I'sJ-iJ's, mi2, m\^, Mj, M'j) = 



j dp3'(f{pi2',P3',a' 



Psq sin e„ 



b3-§l 
0q+arccos Z3 

[ sin 03(193- — Sms'msa 

J sm 0)3 sm U3 I 

1 1 

(/3,Mj -mi2 - ms3, -,ms3|/3, -, j3,Mj - mia) 



q 

+ [^ll,M'j-m{2-ms'^^P3 + 2)^/3, Mj-mi2-ms3 (P3)]</,3=_<^'3) 



1 3 

= y dp3'(fi{pu\p3\a')X3{p3,p'r^;l3j'^j3,j'^,mi2,m[2,Mj,Mj). (43) 
b3-fl 

with X3{p3,p'^;l3,l'^,j3,j'^,mi2,m[2,Mj,Mj) as a one-dimensional integral over 6*3. The 
kronecker-delta Sjns'^,ms3 comes from the fact that our operator does not act on the spin 
of the third nucleon. Thus, the final expression for the matrix element looks like- 

MiM'j,Mj) = 3 

2 ."^'12 '™12 .*'l2 '^12 ''12 ''12 .is 'is .^3 ,'3 12 

j Pl2''^dpi2' J PudPu J Pldp3V{Pl2,P3,Ot) 

(Pl2 , P12 ' ; /12 , ^'12 , S 12 , S'i2 , jl2 , 2 ) "^12 , "^'l2 ) 12 ) 
P3+f 

j C^Pa V(Pl2',P3',tt') 

b3-§l 

(P3 , P3 ; ^3 , ^3 , h Js^mu, m'i2 ,Mj,Mj). (44) 



2. Further Reduction of One-Body Matrix Elements 

At this point, we would like to recall that the operator has a one-body and a two-body 
part and the kinematics make it further evident that the one-body calculations should be 
simpler and involve fewer integrals. The one-body matrix element is- 

= 3(^/1^(0^^(1) + 0^^(2))|v[/.) = ^((^/|0^^(l,2)|vl/,)) (45) 
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where 0^^(1,2) is the one-body operator in the two-nucleon spin-isospin space. This means 
that only one of the three nucleons is involved in the scattering process and we are gifted with 
yet another three-momentum conserving delta-function. For the purpose of this discussion, 
considering nucleon '1' (see Fig. [7]) as the 'struck' nucleon, we can use- 

^ I Pz' k' _ ^ Ps k 

^Pi2' = Pi2 - |, (46) 

after using Eq. (!33|) . We have to be careful here because in our calculations we use the basis 
of two-nucleon spin and isospin. However, we use the fact that the spin and isospin sums 
are performed only over allowed partial waves and hence we can assume that nucleon '1' is 
the struck nucleon and use Eq. (H^ to relate P12' = pu — |- This is very useful because 
we can now reduce integral X2 in the same manner as X3. After going through all the steps 
described in Sec. Ill B II we arrive at- 

il2 >il2 i'"'l2 i"T-12 ,s'i2 i«12 ,^^2 ''12 lis lis ,'3 ,^3 ,mti2 



Pudpu J P3dp3(p{pi2,P3,a) 

+ f Pl2 + f 

j dps' j dpi2'(p{pi2',P3',a') 



|P3-§| bl2-f| 



where. 



^2{Pl2,Pl2 V12, ^12' ■512, S'l2? il2, Jl2? "^12, "^'12? ^^12) 

^3(^3,^3; ^3, ^3, J3, Js, mi2, m',2, Mj, M'j). (47) 



^2 {Pl2 , Pl2 Vl2 , ^'12 , ■S12 , , jl2 , J12 , mi2 , , mti2 ) 

Sq+arccos 212 

sm 6'i2tt6'i2 



Pi2q sin Oq 

arccos Z12I 

— E 

sin(/)i2 sin6'i2 , 

jll I II I \ll I -I I \ 

l'l2, ^^12 ~ ^^^512, -512, '^■5l2rl2, ^512,^12, '^12J 
(/l2,mi2 - mSi2, Sl2,"^Sl2|^12, S12, jl2,"^12) 

{\^V^^,m'-^^-ms'^^{Pl2 - |)'^ii2,mi2-msi2(Pl2)]0i2=0'i2 
^^l>^^,ra>^^-ms\^^PV2 — |)'^«i2,mi2-msi2 (Pl2)]</,-^2=-</,i2) 

(tl2"^^12KSl2"^^s'l2|C'^'^(l,2)|si2mSi2)|tl2"^tl2) (48) 
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and 



Pl2 + T ~ -^12 
Zl2 = 

Pi2g 

~ zi2 - cos 6^12 cos 

COS 012 = -. (49 

sm (Ji2 sm t^g 

Now that we have the structure of the matrix elements, let us look at the operators. We 
have to translate the one-body operator (from App. and the two-body operators (from 
App. [B]) into the two-nucleon spin-isospin space. 



3. One-Body operator in Two-Nucleon Spin-Isospin Space 

The one-body operator as described in App.[A]can not be used directly. So, using Eq. (lA2p 
we write- 

6-(i, 2) = r« + t(?; = (Af + Aftf ) (50) 

i=1...6 

where ti are the operators involving the photon polarization, momenta and nucleon spin and 
the superscripts '(1)' and '(2)' refer to specific nucleons. Then defining 

A = A^'^ + A^'''r, (51) 

where, the superscripts (IS) and {IV) refer to isoscalar and isovector pieces of Ai respec- 
tively, we can rewrite Eq. fl50|) as- 

o-(i,2) = 

i=1...6 

*+(^"-rf )(«.'■' (52) 

And we have, 

{t[,mt,2\{4'^ + rf )|ti2mti2) = 2 (53) 

only for t'^^ = ^12 = 1 and mti2 = 1. (We do not consider the case when mti2 = — 1 because 
it is not possible for a ^He nucleus.) Also, 

{t'i2^h2\i4'^ - rP)\h2mtu) = 1 (54) 

for (10| |00) and (00| ^ |10) only and the rest are zero. 

For the spin part of the operator we have either {t^p + tf^) or {tl^^ — tl^^). And we can 
write- 

(tf ^ + tf ^) = 2(ti + t2)+2i S-Vi (55) 

j=3...6 

where S is the total spin operator of the two-nucleon system and is a vector that contains 
the photon polarization and momenta. Similarly, 

i=3...6 

This operator induces the transitions (00 1 <— |lmsi2) and (Ims'j^al ^ 



19 



4- The One-Body Matrix Element 



At this point, lets summarize how the different pieces combine to give Ai^^{Mj, Mj). 
Since, l3{p3,p'^;l3,l'^,j3,j'^,mi2,m[2,Mj,Mj) does not have any pieces of the Comp- 
ton operator, it is calculated separately. Then using the information from Eqs. ( 153|1 
and (IMI) we calculate the isospin matrix elements in two-body isospin space. Now, 
for each of the non-zero isospin transitions, the spin transition matrix elements are 
calculated using the information from Eqs. (155!) and (156!) . i.e. by calculating the 
expectation values of S or {a^^^ — a^^^) between two-nucleon spin states as dic- 
tated by Eq. (13^ . Now that we have {t[2iT7ti2\{s[2'rns[2\0^^{l,2)\si2msi2)\ti2rnti2), 
we put this in Eq. fHSjl and numerically calculate 22{Pi2,Pi2'', h2, ^12^^12, s[2,ji2,ji2y 
17112, m'i2,T7it 12). Then we plug in the values of X3{p3,p'^]l3,l'^,j3,j'^,mi2,m'i2,Mj,Mj) and 
2^2(pi2,Pi2 Vi2,^i2?'5i2,s'i2, ji2, ji25"^i2,"^i2?'"^^i2) in Eq. (jU]), evaluate the integrals over 
Ps') PS) P12' and pu after introducing the wavefunctions projected into the equivalent ba- 
sis and finally, sum over all angular momentum and isospin quantum numbers to obtain 
M^^{M'j,Mj). 

Next, we shall describe the procedure for the two-body operators and put everything 
back together to obtain the two-body matrix element. 



5. Structure of Two-Body Matrix Elements 
The two-body matrix element can be written as- 
M'''{M'j,Mj) = 3 Yl 



Pl2dpi2 j Pu'^dpu' J pldp3'f{pi2,P3,a) 
dp3 V(Pi2',P3',a') 



I2{P12,P12 V12, l'l2^ -512, S[2,jl2j'l2^ ^12, ^2. "^^12) 

^"3(^3,^3; ^3, ^3, J3, Js. mi2, m[2, Mj,Mj). (57) 
where, l2{Pi2,Pi2 '] ^12, ^'i2> ^u, s[2, ju, j'12, "^12, "^12, "^^12) in its entire glory is given by- 

I2{P12,P12 /12, I'u^ S12, S[2,jl2,j[2^ ^12, m'^2, mtu) = 

dpi2' j Pl2y!'^^,s',,j[,,rn'JPl2') 
{pi2']t[2mti2\d{l,2)\pi2-,ti2mti2)yi,,,si2dl2,mi2{Pl2) 
= j dpi2 ' j Pl-2 Y (^12 

, Tni2 mS 12, Si2, '^■5]^2 1^12 5 '^12? il2? ^12) 

msi2 ,m'^2 

(/l2,mi2 - mSi2, Si2,mSi2|/l2, Si2,jl2,mi2) 
Yll2,m[2-ms[2 (^12 ')^l2,mi2-msi2 (^12) 

(pi2 t[2mti2s[2ms[2\0^^ {1, 2)|pi2; ti2mti2Si2msi2) (58) 
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Let us then see how (pi2 ^ ^I2"^^'i2'^i2"^'^i2l^^ (1' 2)|pi2; ti2"^^i2'5i2"^'5i2) is evaluated. 
6. Two-Body Operator in Two-Nucleon Spin-Isospin Space 

The two-body operator is given by Eq. fIBll) and each of the terms in that equation are 
expressed in the equations following it. Symbolically, we can write each of those diagrams 
as (for example, for the first diagram in Fig. [S])- 

(a) _ _^9a /-.(I) . -.(2) _ J1)^(2)n ^-^^ ^' 



■^-fNN — f2 ^ '3 '3 



= F^^Kf^^^ ■ r(2) - ri^Vf )e ■ a' t' ■ a' (59) 
In general, the spin-isospin structure of the two-body operator would be- 

The only transitions in isospin space that are non-zero are (00 1 ^ |00) and (10 1 ^ |10) and 
both of these are- 

(00|(f « • r - ri^Vf )|00) = -(10|(r « • f - rj'Vf )|10) = -2. (61) 
For the spin part of the operator we can write- 

a'-Aa'-B = ]^{ds + da) (62) 



where Os is the symmetric combination of the LHS of Eq. fl62|) and Oa is the asymmetric 
combination under the interchange of nucleons. Then, we can write Os as- 

ds = ■ AS ■ B -2A- B -2iS ■ {Ax B). (63) 

Here, 25" = + c?^ and we have used the identity a ■ A a ■ B = A ■ B + ia ■ {A x B). 
Since this is a spin-symmetric combination Og induces the transitions (00 1 ^ |00) and 
(Ims'^jl ^ 11^-512). The different spin transition matrix elements can be written in terms 
of general vectors A and i? as - 

(00|a|00) = -2A-B 
{Ims'^^lO slims 12) = -2A- BSms[2,msi2 

+8 ^(1, a, 1, 1, 1, ms[^)il, p, 1, msi2|l, 1, 1, Ai)^^ 5^ 

fiaf3 

+i2V2^(l,7,l,msi2|l,l,l,ms;2)(lx 5)t (64) 
7 

where we have introduced spherical components of vectors A and B. Similarly, we can 
simplify the asymmetric combination to- 

da = 2S- A{a' -a^)-B- i{a' - a^) -{AxB). (65) 
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This, being a spin-asymmetric combination, it induces transitions between (00| ^ |lmsi2) 
or {lms'i2\ ^ |00). After going through some spin algebra we obtain- 

(Wi2|a|00) = 2v^^(l,a,l,/5|l,l,l,m^)At5;f 

a/3 

-tY,^ms[„yi^y<B)\. (66) 

7 

We can calculate {00\Oa\^Tnsi2) by taking the hermitian conjugate of Eq. (l66i) . 
7. The Two-Body Matrix Element 

We can now summarize how the different pieces combine to give Ai'^^{Mj, Mj). Again, 

^3iP3yP3jhJ3,j3,j'3,nT'i2,m'i2y^Jy^j) ^ocs uot havc any pieces of the Compton opera- 
tor, it is calculated separately. The isospin transition matrix elements are simple in this 
case and given by Eq. (l6T|l . For each of these two isospin transitions, the spin transi- 
tion matrix elements are calculated using the expressions from Eqs. fl64l) and fl66p . (Re- 
member, the vectors A and B are different combinations of photon polarization and mo- 
mentum vectors and nucleon momenta.) Now that we have (pi2 ^'i2'^^i2Si2^'S']^2 

10^^(1,2) 

\'Pi2',ti2mti2Si2rnsi2), we put this in Eq. ( !58|) and numerically calculate it. Then we plug in 
the values of Is{p3,p'3;h,l'^, js, j'^,mi2,m[2, Mj, M'j) and l2iPi2,Pi2';h2,l'i2^ s^, 3^2, ji2j'i2^ 
mi2, '"^^12) in Eq. ( 1571) . numerically evaluate the integrals over p^', p^, pi2 ' and pi2 after 
introducing the wavefunctions projected into the equivalent basis and finally, sum over all 
angular momentum and isospin quantum numbers to obtain Ai^^{Mj, Mj) (Eq. (1571) ). 
Finally, we add the one-body and the two-body pieces to obtain- 

M{M'j, Mj) = M'^iM'j, Mj) + M^^{M'j, Mj) (67) 

and using this we can now calculate various observables. In the next section we discuss the 
observables that we concentrate on and then in Sec. |V] we discuss some of our results for 
those observables. 

III. OBSERVABLES 

A. Differential Cross-section 

The differential cross-section (dcs) is proportional to the spin-averaged sum of the square 
of the scattering amplitude: 

Here, Mj and M'j are the initial and final spin projections of the nucleus respectively, M^He 
is the mass of the ^He nucleus and s{ijj) is Mandelstam s. The factor of \ comes from the 
averaging over the initial ^He spin and photon polarization states. 



22 




FIG. 10: The double-polarization observables. The inset in the bottom left shows possible circular 
polarization states for the incoming photons (on the left) and possible initial target spin polarized 
states (on the right). 



da_\ / dcr \ 

= . )t:' - (69) 



B. Double-Polarization Asymmetry 

The double-polarization asymmetry involves circularly polarized photons and a spin- 
polarized target (see Fig. [TU]) . The expectation is that these observables can provide insight 
into the spin polarizabilities of the target. 

When the target is polarized along the beam direction (parallel or anti-parallel), the 
corresponding observable is called the parallel target polarization asymmetry and is defined 
to be: 

' da \ I i dcr 

Parallel (anti-parallel) arrows in the subscript symbolize target polarization parallel (anti- 
parallel) to the beam helicity and A denotes the helicity of the incoming photon. This 
observable can be depicted only in terms of the numerator of Eg. (1691) as a difference in 
cross-section as- 

The target may be polarized along ±x too. In this case, the observable is called the 
perpendicular polarization asymmetry and is given by: 

I dcr\ ( — \ 

^-,(A=±i) - Td^-—7d^- (71) 

Again, the direction of the second arrow in the subscript denotes the target polarization is 
along the +x or —x directions. A is the helicity of the incoming photon. For this observable 
one must be careful in defining the matrix elements because the spin state of the target is 
an eigenstate of Sx and not of Sz- The eigenstates of Sx should be first expressed in terms 
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of those of Sz and then one can proceed with the calculations. As before, this observable 
can also be expressed as a difference in cross-sections as- 

It should be mentioned here that the results reported in this paper are for and A^;. 



IV. HOW DO POLARIZABILITIES ENTER INTO THE CALCULATION? 

The one-body amplitude (App. |A]) has six invariant structures Ai . . . Aq that depend on 
the photon energy. When these amplitudes are Taylor-expanded in u around ~ we 
obtain the following expressions [l^ 

As = i^[Z{Z + 2h) -{Z + COS0] + Af + (71 - (72 + 274) cos{e))u^ + ^(o;^), 



Z + kYuo 



A, = -^^^ + Af + 73-^ + 0{.% (73) 

Here, Af , Af and Af are contributions from the 7r°-pole graph (see Fig. [31 lower-most 
diagram). The zeroth-order term in uj is the leading-order Thomson term (as in Eq. ([7])) 
and the first-order term contains contributions from the anomalous magnetic moment. The 
uj"^ term depends on a and (3 and the spin polarizabilities (7's), appear in the term. At 



0{e'^Q) the following results for the polarizabilities can be obtained j59|, l60| by matching 
Eq. dMl) to Eq. ([73]): 



12.2 X 10"^ fm^ 

2„2 



3847r2/>. 



Pp = Pn= ^J^f, = 1.2 X 10-^fm^ 



2 2 

^ip = = 0/3/2 2 = 4-4 X 10-Hm\ 

72p = l2n = ^J'ff, , = 2.2 X 10-Hm\ 

2 2 

73p = 73n = , = 1.1 X 10-^fm^ 

74P = 74n = - .jYf, , = -1.1 X 10-^fml (74) 
3847r'^ j^m^ 
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The manner in which we investigate the impact of the neutron polarizabihties is by 
varying them when calculating any observable around these central 0{e^Q) values. We 
introduce six new parameters as corrections to the 0{e'^Q) values, which we call Aa„, A/3„ 
and A'y,r^{i = 1,2,3,4): 

A2 = A^iOiQ')) + APnCO^ 

As = As{0{Q')) + (A7i„ - (A72„ + 2A-f,^) cos{9))u;^ 

= A4(c»(g3)) + A72„^^ 

A, = A5(0(g')) + A74„cu3, 

A, = A(0(g3)) + A73„cu3. (75) 

But, note that, this only includes one set of higher-order contributions. It does not represent 
a full calculation at 0{e'^Q^) and it only examines a limited set of the contributions generated 
by the A-isobar. 

For one particular plot only one of these parameters, for instance, Aa„ is varied with the 
rest being set arbitrarily to zero. This gives us a measure of the sensitivity of that particular 
observable to Aa„. In principle, one should fit these parameters to experimental data and 
thus extract them, but at present we are still waiting for data on most of these observables. 

V. RESULTS 

A. Coherent Compton Scattering 

In Fig. [TT]we plot the dcs for coherent 7^He scattering. For the results reported in this 
section we employ the ^He wavefunction obtained using the Idaho chiral N^LO NNpoten- 
tial [so'] with the cut-off at 500 MeV together with a chiral three-nucleon interaction [52I]. All 
the different panels are for calculations at different energies (60, 80, 100 and 120 MeV) in the 
cm. frame. Each of the panels shows the dcs calculation at different orders- 0(e^), Impulse 
Approximation (lA) and 0{e^Q). Here, Impulse Approximation means that the calculation 
is done up to 0{e^Q) but does not have any two-body contribution- this is not the full 
0{e'^Q) calculation. As expected, we see that there is a sizeable difference between the lA 
and the 0{e^Q) dcs which means that the two-body currents are important and cannot be 
neglected. Also, we see that the difference between the O(e^) and 0{e^Q) is very small at 
60 MeV and it gradually increases with energy. This can be attributed to the fact that, as 
energy increases, the energy-dependent contributions from the Ai . . . A^ terms increase and 
this means that there is an opportunity to extract the neutron polarizabihties. At the lowest 
energy the dcs is still dominated by the proton Thomson term and hence the 0{e^Q) dcs 
is closest to the 0{e^) one, or in other words, xPT converges well there. Another notable 
feature is that the 0{e^) dcs is independent of the cm. energy at forward angles, or we can 
say that the dcs at 0{e^) converges to a single value at zero momentum transfer. 

Now, focusing just on the O^e^Q) curves for all four energies it is evident that there is a 
rapid decrease in the dcs as one goes from 60 MeV to 120 MeV. It remains to be seen whether 
this drop-off is compensated by addition of the A-isobar to the theory. Also notable is the 
fact that the ^He Compton scattering dcs is around four times bigger in magnitude compared 
to deuteron Compton scattering dcs. We conclude that Compton scattering experiments on 
^He are potentially more accurate than on deuteron. 
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FIG. 11: Comparison of differential cross-section calculated at different orders in the cm. frame. 
The top two panels are for calculations at 60 MeV (left) and 80 MeV (right) and the bottom two 
panels are calculations at 100 MeV (left) and 120 MeV (right). The dashed (blue) curves are the 
0{e'^) results, the solid (black) curves are the impulse approximation results and the dot-dashed 
(red) curves are the 0{e^Q) results. 



Next in Figs. [I2] and [I3] we plot the 0{e^Q) dcs at different energies by varying Aa„ 
and Af3n respectively. Here we vary Aq;„ between (—4. . .4) x 10~^ fm^ and A/3„ between 
(—2 ... 6) X 10~^ fm^. (The range of variation of Aj3n is chosen in this manner because of 
a strong paramagnetic contribution from the A-isobar jssl . |37|.) What is striking in these 
figures is that because the ^He dcs is large compared to that of the deuteron, the sensitivity 
to Aan and A/3„ is also similarly magnified. It should be evident from the figures that the 
size of the absolute differences in the dcs as we vary Aa„ and A/?„ is roughly the same at 
all energies. However, we would like to advocate that if such measurements were necessary 
(we expect that a„ and /3„ could be determined from ongoing coherent 'yd experiments at 
MAXLab at Lund) then it would be better to do them at lowest energy possible. At the 
lower energies, the dominant contributions come from a„ and apart from the proton 
Thomson terms. In this regime, we best understand the dynamics and the extraction will 
not be tainted by contributions from the spin polarizabilities entering via A3 . . . Aq and other 
higher-order contributions. The figures show that even at 60 MeV, the sensitivity to a„ and 
Pn is sizeable. Another message to take away from these two figures is the beautiful feature 
in Fig. [13] that sensitivity to /3n vanishes at = 90 deg and the curves themselves turn over. 
Thus, in principle one could extract a„ and /5„ independently from the same experiment as 
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FIG. 12: The differential cross-sections in the cm. frame with varying Aa„. The top two panels 
are for calculations at 60 MeV (left) and 80 MeV (right) and the bottom two panels are calculations 
at 100 MeV (left) and 120 MeV (right). The solid (black) curves correspond to the fuh 0{e'^Q) 
results. The long-dashed (blue) curves correspond to Aa„ = —4 x lO^^fm'^, dot-dashed (red) to 
Aa„ = — 2 X 10~"^fm^, dotted (magenta) to Aon = 2 x lO^^fm'^ and dashed (g reen) to Ao^ — 
4 X W~Hm^ 



opposed to measuring one and then using the sum rule to get the other. A measurement 
near = 90 deg would yield a„ and then one could perform a forward/backward dcs ratio 
measurement to extract The various curves for Aa„ suggest that this ratio should not 
vary much as Aa„ is varied, but, on the other hand the effect of A/5„ would be magnified 
compared to what we show in Fig. [131 

In summary, the extraction of and /3„ is possible through coherent Compton scattering 
on unpolarized ^He. Extraction of a;„ and j3n in this manner would also serve as a check for 
the extractions of a„ and /5„ from the unpolarized 7d experiments at MAXLab. 



B. Polarized Compton Scattering 

In Fig. [H] we plot A^ vs. the cm. angle at 120 MeV and the different panels correspond 
to varying the four spin polarizabilities one by one. Here we choose to vary 7j„, i = 1 . . .4 
between ±100% of their C(e^Q) prediction (Eqs. fl7^ ). Looking at the plots, the message 
one gets is that this observable is quite sensitive to 7i„, 72,1 and 74„. Especially with the 



27 




FIG. 13: The differential cross-sections in the cm. frame with varying A/3„. The top two panels 
are for calculations at 60 MeV (left) and 80 MeV (right) and the bottom two panels are calculations 
at 100 MeV (left) and 120 MeV (right). The solid (black) curves correspond to the fuh 0{e^Q) 
results. The long-dashed (blue) curves correspond to A/3„ = — 2 x lO^^fm'^, dot-dashed (red) to 
A/3„ = 2 X 10~'^fm^, dotted (magenta) to A/3„ = 4 x 10~'^fm^ and dashed (green) to A/3„ = 
6 X lO-'^fm^ 



expected increase in photon flux at HI7S such differences in cross-sections should easily be 
measured and we can expect to extract a linear combination of 7i„, 72^ and 74^. The linear 
combination that readily comes to mind is 70^ or 7^^ but we should bear in mind that we 
are measuring as a function of angle which means that we should be able to extract the 
combination 7i„ — (72^ + 274„) cos6' (see Eq. (!73i) ). Observing the different plots in Fig. [H] 
we can already see the effect of the cos^ term in the panels for 72^ and 74^. 

Next, in Fig. [15] we plot A^; vs. the cm. angle at 120 MeV and the different panels 
correspond to varying the four spin polarizabilities one by one. This figure also suggests 
that we are sensitive to a combination of the same spin polarizabilities as in but this 
combination is obviously not the same as before i.e. 7i„ — (72,1 + 274„)cos6'. We can say 
that because the curves in the right panels of Fig. [15] should coincide at 90 deg because 
of the cos 6' term but they do not. Hence, we should be able to extract a different linear 
combination of the 7„s from A^,. 

Thus, we are sensitive to two different linear combinations of 7i„, 72^ and 74^ through A^ 
and A^: and at least we can expect that we can extract one, say 71^ unambiguously, out of 
the three and put constraints on the values of the other two. For scattering we found that 
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FIG. 14: The four panels above correspond to with varying each of A7i„ (top left), A72n (top 
right), A73„ (bottom left) and A74„ (bottom right), one at a time. The calculations are done in cm. 
frame at 120 MeV. The solid (black) curves correspond to the full 0{e'^Q) results. The long-dashed 
(blue) curves correspond to A'jin = — 7m(0(e^(5)), dot-dashed (red) to A7j„ = — 7j„(0(e^(5))/2, 
dotted (magenta) to A.'jin = 7i„(0(e^(5))/2 and dashed (green) to A7j„ = 7m(0(e^<5))- 



we were sensitive to a combination of 7i„ and 73^ 3^ . This means that if we combine all the 
information, we should be able to extract at least two of the four spin polarizabilities and 
constrain the remaining two. Also we can definitely hope that through polarized Compton 
scattering on deuteron and ^He we can perform a nuclei independent extraction of at least 
one of the neutron spin polarizabilities, presumably 7i„. 

It is worth mentioning here that the curves in Figs. [T3] and [T5] look very similar to those 
obtained if one calculates the reaction jn —>■ jn, i.e. Compton scattering off a free neutron. 
This would then suggest that polarized ^He indeed behaves as an "effective" neutron target. 



VI. POLARIZED ^HE IS INTERESTING 

In this section we would like to point out a couple of very interesting facts about Compton 
scattering on polarized ^He. 

1. We know that ^He is a spin-| target and hence, we should be able to write the photon 
scattering amplitude as a sum of six structure functions like in Eq. (IA2|) . However in 
this case these six functions have to be a sum of one-body and two-body parts and 
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FIG. 15: he four panels above correspond to A^; with varying each of Ajin (top left), A72n (top 
right), Aj^n (bottom left) and A74„ (bottom right), one at a time. The calculations are done in cm. 
frame at 120 MeV. The solid (black) curves correspond to the full 0{e'^Q) results. The long-dashed 
(blue) curves correspond to A'jin = — 7j„(0(e^(5)), dot-dashed (red) to A7j„ = — 7j„(0(e^(5))/2, 
dotted (magenta) to A'jin = 7i„(0(e^(5))/2 and dashed (green) to A7j„ = 7m(0(e^<5))- 



are not as simple as the AiS in Eq. (IA2I) . 

i=i...e 

Al"^ = Al^ + Af. (76) 

Here, A}^ are the same as in Eq. flA2p to the extent that ^He is an "effective" neu- 
tron. Meanwhile, we found that, numerically, Aj^,i = 3... 6 are negligible. Hence 
we conclude that Aj^ does not contribute to the spin structure functions at 0{e'^Q). 
Moreover, since the two protons (to a good approximation) have spins anti-aligned in a 
spin-polarized ^He target, we can safely assume that all the sensitivity to the spin po- 
larizabilities in the observables come from the unpaired neutron alone. This reasoning 
supports our claim (made in the previous section) that the curves for ^He Comp- 
ton scattering resemble those obtained from n-Compton scattering. Thus, Compton 
scattering on '^He is an ideal avenue to extract the neutron spin polarizabilities. 

2. Let us consider the double polarization observables, and A^;. If we take Eq. fl76|) 
to calculate these observables then we shall obtain expressions similar to Eqs. (3.19) 
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and (3.22) of Ref. 13] but with Ai replaced by A^^^. From these equations we know 
that to the lowest order in u {0{uj^)) the spin polarizabilities manifest themselves via 
the interference terms of . . . Aq with Ai. Hence, to the extent that polarized ^He 
is an effective neutron, we can expect the A3 . . . Aq for the neutron will be multiplied 
by Ai^'^, which- at least at low energies- is dominated by the two proton Thomson 
terms, thus giving a much enhanced effect compared to 7d. 



VII. SOURCES OF UNCERTAINTIES 
A. Wavefunction Dependence 

In calculating the matrix elements (Eq. ( l24|) ). ideally, the wavefunctions and the operator 
should be calculated within the same chiral framework, so that the formalism is consistent. 
Although few-nucleon bound state wavefunctions based on chiral nuclear interactions have 
been calculated, the currently used approaches are not entirely consistent with the approach 
used for the Compton scattering operators. Therefore, we employ several wavefunctions 
based on chiral perturbation theory and additional wavefunction obtained from the state- 
of-the-art interaction models. In this way we expect to cover the entire dependence on the 
choice for the wavefunctions. Specifically we chose- 



1. The phenomenological potential Al^lS jSJJ with the 3A^ interaction Urbana IX [53 



2. The Idaho N^LO chiral potential [50|] with the cut-off A = 500 MeV with a chiral 3A^ 



force j52| (called 3NFB here). 



3. The Idaho N^LO chiral potential [SOl] with the cut-off A = 500 MeV with another 
chiral 3A^ force [s^] (called 3NFA here). The difference between 3NFA and 3NFB 
is that they are parameterized differently so that they reproduce the triton and alpha- 



particle binding energy but lead to different spectra for heavier nuclei |52l ]. 

4. The Idaho N^LO chiral potential Q with the cut-off A = 500 MeV and no 3NF. 

5. The Idaho N^LO chiral potential [soj with the cut-off A = 600 MeV and no 3A^F. 



6. The NLO chiral potential with dimensional regularization |5J] with the cut-off in the 
Lippman-Schwinger set to A = 600 MeV and no 3NF. 

In the future, it will also be interesting to employ wavefunctions based on NLO, N^LO 



and N^LO interactions of 55||, which use a different regularization scheme and allow one 
to compare several orders of the chiral interaction within the one framework. However, 
we believe that the set of interaction from above should give a reasonable idea of possible 
wavefunction dependence in this exploratory study. What we want to demonstrate through 
the choice of these particular wavefunctions is that- 

• A contrast between the NLO chiral wavefunction (A = 600 MeV) and the Idaho 
N^LO wavefunction (A = 600 MeV) will demonstrate the effect of terms of higher 
chiral order. 

• A contrast between the Idaho N^LO wavefunction (A = 600 MeV) without a 3NF and 
the Idaho N^LO wavefunction (A = 500 MeV) will demonstrate the effect of varying 
the cut-off which translates into studying the impact of short-distance physics. 
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• A contrast between the Idaho N^LO wavefunction (A = 500 MeV) and the Idaho 
N^LO wavefunction (A = 500 MeV) with ?,NFA will demonstrate the effect of 'a' 
three-nucleon force. 

• A contrast between the Idaho N^LO wavefunction (A = 500 MeV) with "iNFA and 
the Idaho N^LO wavefunction (A = 500 MeV) with 3NFB will demonstrate the effect 
of choosing three-nucleon forces with different parameters. 

• A contrast between the Idaho N^LO wavefunction (A = 500 MeV) with 3NFB and the 
AV^18+Urbana-IX wavefunction will demonstrate the effect of choosing a wavefunction 
constructed with a phenomenological potential. 

We have calculated the dcs at two different energies (60 & 120 MeV) and also A^. and 
at 120 MeV with these six wavefunctions. The results are plotted in Fig. UM 




FIG. 16: The top two panel show the dcs at 60 MeV (left) and 120 MeV (right) calculated with 
different wavefunctions. The bottom left panel corresponds to at 120 MeV and the bottom right 
panel corresponds to A^ at 120 MeV. These have also been calculated using different wavefunctions. 
The dashed-shortdashed (grey) curves correspond to j4V18 + UrblX, solid (black) corresponds to 
Idaho N^LO + 3NFB, the long-dashed (blue) curves correspond to Idaho N^LO + 3NFA, the 
dot-dashed (red) curves correspond to Idaho N^LO with A = 500, the dotted (green) curves 
correspond to Idaho N'^LO with A = 600 and the dashed (magenta) curves correspond to NLO 
chiral wavefunction with A = 600. All the calculations have been done in the cm. frame. 

Let us first compare the top two panels that show the dcs. The first thing to notice 
is that the shape of the dcs is the same no matter what wavefunction we use. Another 
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fact that stands out is that all the chiral potentials produce dcses that are close to one 
another, whereas the choice of AV^lS+Urbana-IX causes the dcs to be quite different. Also 
the size of the difference in the dcs due to the choice of different wavefunctions increases as 
we go from 60 to 120 MeV. One major difference between the chiral wavefunctions and the 
y4V^18+Urbana-IX wavefunction is that the latter has explicit high-momentum components 
in the NN interaction and this causes the dcs also to be significantly different. Having said 
that, it is essential to acknowledge the fact that the dcs is sensitive to the choice of different 
chiral potentials too. Of these, the chiral wavefunctions with three-nucleon forces and the 
N^LO chiral wavefunctions are actually of a higher chiral order than our Compton scattering 
operator and we expect some sensitivity to the higher chiral order terms. Nevertheless, 
looking at the dcs it can be said that the uncertainty due to the choice of wavefunctions 
(chiral wavefunctions only) is < 15% at 120 MeV. Of course this is lower at 60 MeV (< 8%) 
and we advocate that if a measurement of the dcs is necessary to extract a„ or /5„ it should 
be done at lower energies. 

Next, observe the lower panels of Fig. [121 What immediately stands out is that even 
though these calculations were done at 120 MeV, the sensitivity to the choice of wavefunc- 
tions is much smaller than in the dcs. This is a good sign because we want to use these 
double-polarization observables to extract the 7„s. The uncertainty due to the choice of dif- 
ferent wavefunctions in these observables is < 7.5%. Also, the calculations with the different 
chiral potentials converge to a single value {i.e. uncertainty is ~ 0%) at forward angles for 

and this augurs well for us because (as we see in Fig. [T^ there is maximum sensitivity 
to the 7„s at forward angles. Similarly, a measurement of at forward angles will also 
reduce the effect of the wavefunction dependence. But, it should be pointed out that, ex- 
perimentally, at forward angles there are large background effects and it will be difficult to 
separate the 7^He scattering events from the background at HI7S. 



B. Position of Pion-production Threshold 

In HB^PT it is assumed that the nucleon rest four momenta dominate their total four 
momenta, in other words, the nucleons are virtually at rest. In consequence the pion- 
production threshold occurs at a; = m^r, because the nucleus cannot recoil. Hence, this 
assumption creates a problem for all processes involving baryons near the pion-production 
threshold. The position of the pion-production threshold in 7 ^He scattering is actually at 
a photon energy of 

2 

uj = -B.He + m^ + — ^ + 0{ml), (77) 

llVhHe 

and not at = m^r. Hence, for the purpose of understanding the effect of this in our 
calculations, we redefined 

as the energy 'going into' the 7A^ amplitude and calculated the dcs at uj =100 MeV. Bzue is 
the binding energy of the '^He nucleus in Eq. (1781) . We found that there was a ~ 5% difference 
in the dcs at 100 MeV. In A^,. the difference was ~ 1% and in A^ the difference was ~ 3% 
at 100 MeV. It is encouraging that the percentage uncertainty in the double polarization 
observables is again less than that in the dcs. This investigation was just exploratory and a 
proper treatment of this issue is an important topic for future study. 
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C. Boost Corrections 



Our calculation of the 7'^He scattering process is in the 7^He cm. frame. However, 
the one-body amplitude (Eq. flA2p ) is defined to be in the 7N cm. frame and similarly the 
two-body amplitude is defined to be in the 7NN cm. frame. Thus, we need to boost these 
amplitudes to the 7NNN cm. frame. This is necessary to ensure that the vectors e and e' 
stay orthogonal to k and k' respectively in the new frame. Consequently, a 'boosting' term 
has to be added to the 7N amplitude which is given by- 



O 



(IB) 
Boost 



1 



^e' ■ke-k' + -(e' ■ ps e ■ k' + e ■ ps e' ■ k) 



1 



2M2cj 



(1) 



(79) 



Here, all symbols have their usual meanings. The structure of the two-body operator is 
such that it is Galilean boost invariant. We calculated the effect of this boost at 100 MeV 
and found that it is truly a perturbative effect and the difference is < 0.01% in the dcs. 
This is because the boost corrections to the transition matrix elements that dominate are 
really very small if not zero. The boost pieces contribute to those matrix elements that are 
themselves small and hence play only a small role in the observables. 



VIII. SUMMARY AND FUTURE PROSPECTS 

We have performed the first 0{e^Q) HB^PT calculations of Compton scattering on ^He. 
Indeed we are not aware of any previous theoretical calculation of this reaction process in 
any framework. We have shown that ^He Compton scattering is an exciting prospect for 
mining the neutron polarizabilities. Polarized ^He is unique and interesting because it 
does seem to behave as an "effective neutron" (refer to Sec. IVII) in Compton scattering. 
With the added advantage of a significant Compton cross-section due to the presence of the 
two protons, whose Thomson terms can interfere with the neutron polarizabilities, we get 
a significantly magnified effect compared to 7d. The results reported in Sec. |V] are quite 
encouraging because they show that not only can we access a„ and /?„ through the dcs but 
we can also hope to extract two different linear combinations of 7i„,, 72™ and 74^ from the 
double-polarization observables, and A^,.. 

Meanwhile, priorities for HI7S (at Triangle Universities Nuclear laboratory) were dis- 
cussed at a Compton workshop held prior to the Fifth International Workshop on Chiral 
Dynamics, 2006 [56] and they included polarized 7p, "yd and j^He experiments. The idea is 
that the 7p —>■ jp experiments will reveal 7ip . . . 74^ and the set for the proton polarizabil- 
ities will be complete. This will set the stage for the neutron polarizabilities. Moreover, if 
the MAXLab experiments generate enough data to extract a„ and Pn unambiguously, then 
7^He can be used to focus on the neutron spin polarizabilities. One can also hope that 
coherent ^He Compton scattering measurements may be done at MAXLab in the future. 
The expectation is that, at the completion of all these planned experiments we will be able 
to complete the set of polarizabilities for the neutron. 
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However, there is a lot of scope for improving the theoretical calculations on ^He Comp- 
ton scattering. This being the first calculation was exploratory in nature. Improvements 
on the theoretical side could include higher-order (O(e^Q^)) calculations and/or explicit 
inclusion of the A— isobar. (We have seen that explicit inclusion of the A-isobar causes no- 
ticeable changes in the results for '^d scattering.) Also dedicated effort is required to shrink 
the theoretical uncertainties associated with the choice of nuclear wavef unctions. Hilde- 



brandt etal. |36l. 1371] have demonstrated that by resumming the intermediate NN scattering 



states they were able to restore the Thomson limit for deuteron Compton scattering and 
also significantly reduce the dependence on the choice of deuteron wavefunctions. Such an 
exercise can also be done for ^He Compton scattering. Effort is also required to rigorously 
study the effect of the incorrect position of the pion-production threshold due to the "static 
approximation" for nucleons in HB^PT . These effects have been studied for ird scattering 
by Baru etal. 57] and for 76? — > Tc^nn by Lensky etal. [H^]- Using similar principles, we 
can endeavor to resolve this issue for Compton scattering on nuclear targets. 

To conclude, it should be reiterated that we have taken the crucial first step. Of course, 
much needs to be done but we are confident that with sustained efforts we shall be able to 
extract the neutron polarizabilities from forthcoming 7^He scattering data. 
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APPENDIX A: 7N AMPLITUDE 



After evaluating all the diagrams in Figs. [3] and HI it becomes evident that the photon 
scattering amplitude T^tv is a sum of six operator structures that can be expressed as: 



T^^ = e^{Aie' ■ e + A2e' -ke-k' + iA:ia ■ {e' e) 

+iAia ■ {k' X k)e' ■ e + iA^a ■ [{e' x k) e ■ k' - {e x k') e' ■ k] 
+iAQa-[{e' X k')e-k' -{ex k)e' -k]} 

6 

= ^ ] Aiti 

i=l 



where, 



(Al) 
(A2) 



ti 
ti 
ti 
ti 

h 

ta 



eH' -ke-k' 
ie^a ■ (e' X e) 



ie^a ■ {k' X k) e' ■ e 

ie^&-[{e' xk)e-k' -{ex k')e' ■ k] 

ie^a ■ [{e' X k') e ■ k' - {e X k) e' ■ k] 



(A3) 



The six invariant amplitudes, Ai...Aq, are functions of the photon energy, uj, and the 
Mandelstam variable t. Defining T = uj/m.,^ and t = — 2T^(1 — cos 6), where 6 is the center- 
of-mass angle between the incoming and outgoing photon momenta, one finds 15, 59, [goI . 



A, 
A2 
^3 



A, 
A, 
A, 



M 87r/2 



1 - Vl - T2 
2-t 



2-t 



- arctan — /i(T, t) 



M2 87r/2m^ {-tf/^ 



UJ 



2M2 



[Z{Z + 2k) - {Z + kY cos 6] + 



-t 

[/i(T,t)-/2(T,t)], 

{2Z - l)gAmjr Tt 



1-t 



T 



arcsin^ T - T + 2T^ sin^ ^/3(T, t) 



[Z + nfu 



2AP +47r4>/'^^''^^' 
{Z + nfuj {2Z - Y 



2M2 



[/5(T,t)-2T2cose/3(T,t)], 



Z{Z + k)uj ^{2Z -l)gAU)'^ T 



2M2 



Sn^f^m^ (1 - t) 
[h{T,t)-2T'h{T,t)], 



(A4) 
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where Z = 1(0) for the proton (neutron) and: 

(l-^)v^ 



/i(T,t) = / (iz arctan 
Jo 

T(T,^ - f\ 2(l-z)v/-t(l-T2z2) 

h{T,t) - dz 4(l_^2^2)_^(l_^)2' 

/3(T,t) = dx dz^ ^-^ '- 



^0 

1 



Tz TzS 



'-^ 2(1-2) Tz 
l4{T,t) = I dx I dz — arcsin-^, 



^0 



^0 



l5{T,t) = I dx I dz arcsin— , (A5) 



with 



5* = Vl-T222-t(l-2)2x(l-x), R = ^/l-t{l-zyx{l-x). (A6) 



APPENDIX B: 7NN AMPLITUDE 

The two-body diagrams that contribute to the Compton scattering process at 0{e^Q) are 
shown in Fig. [5l The two-body amphtude can be expressed (in the 7NN cm. frame) [isl . 
3 113 as: 

32^2 
2/, 

where 



T^2B^ = (f 1 . r 2 - rgVI) (t(») + + + t^'^ + &^), (Bl) 



2[cu2-m2 -(p-p' + i(A; + A;'))2] 



+ (1 ^ 2), (B2) 



_ e ■ ■ (p - - |(fc - k'))a' ■ {p- p' + |(fc -k')) 

2[{p- p' -\{k- k')Y + ml][{p- p' + \{k- k')Y + ml] 
+ (1 ^ 2), (B3) 
e' -{p-p' + \k) a' -ea' ■ [p-p' + l{k-k')) 



lu)'^ — ml 

I TV 


— {p — p' 


+ \ik + k')Y][{p- 


-p' + 


i(fc-fc'))' + <] 


(1 ^ 2), 










e ■ 


(p - p' + 


\k') a'-{p-p'- 


lik- 


-k'))a^-e' 


[ct;2 — m2 


— (p — 


+ l{k+k')mp- 


' p' — 
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(1 ^ 2), 










2e ■ (p — J?' 


+ \k')e' 


■ {p- p' + \k) ■ 


{p- 


p' -\{k-k')) 



(B4) 



(B5) 



[uj^-ml-{p-p' + \{k + k')f][{p-p' -\{k-k')f + ml] 

a'' -ip-p' + Uk-k')) , , 

X ^ ^\ ^+ 1 ^ 2 , B6 

[(p-p' + i(fc-A;'))2 + m2] 
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— * — * 

Here again, e,e',k,k' have their usual meaning. p{p') is the initial (final) momentum of 
the nucleon inside the deuteron, and is twice the spin operator of the first (second) 

nucleon. The numbering of the nucleons is arbitrary and hence (1 2) denotes the term 
when the nucleons are interchanged. 
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